Magnetic anisotropy of epitaxial (Ga,Mn)As on (113)A GaAs 
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The temperature dependence of magnetic anisotropy in (113)A (Ga,Mn)As layers grown by molec- 
ular beam epitaxy is studied by means of superconducting quantum interference device (SQUID) 
magnetometry as well as by ferromagnetic resonance (FMR) and magnetooptical effects. Experi- 
mental results are described considering cubic and two kinds of uniaxial magnetic anisotropy. The 
magnitude of cubic and uniaxial anisotropy constants is found to be proportional to the fourth and 
second power of saturation magnetization, respectively. Similarly to the case of (001) samples, the 
spin reorientation transition from uniaxial anisotropy with the easy axis along the [TlO] direction at 
high temperatures to the biaxial (100) anisotropy at low temperatures is observed around 25 K. The 
determined values of the anisotropy constants have been confirmed by FMR studies. As evidenced 
by investigations of the polar magnetooptical Kerr effect, the particular combination of magnetic 
anisotropics allows the out-of-plane component of magnetization to be reversed by an in-plane mag- 
netic field. Theoretical calculations within the p — d Zener model explain the magnitude of the 
out-of-plane uniaxial anisotropy constant caused by epitaxial strain but do not explain satisfactorily 
the cubic anisotropy constant. At the same time the findings point to the presence of an additional 
uniaxial anisotropy of unknown origin. Similarly to the case of (001) films, this additional anisotropy 
can be explained by assuming the existence of a shear strain. However, in contrast to the (001) 
samples, this additional strain has an out of the (001) plane character. 

PACS numbers: 75.50.Pp, 75.30.Gw, 73.61.Ey 



I. INTRODUCTION 

Since many decades, a lot of attention has been de- 
voted to ferromagnetic semiconductors. More recently, 
the intense research has been triggered by the syn- 
thesis of the (III,Mn)V diluted magnetic semiconduc- 
tor (Ga,Mn)As, 1 which has become the canonical ex- 
ample of a dilute ferromagnetic semiconductor. 2:3 It has 
been demonstrated that a number of pertinent proper- 
ties of this material can be explained by the p — d Zener 
model. 3-6 Magnetic anisotropy of strained (Ga,Mn)As 
layers can be calculated within this theory, and many 
experimental studies 7-15 were devoted to verify its pre- 
dictions. However, despite these intense studies, some 
important features of magnetic anisotropy in this system 
are at present not completely understood. 

An example of such a property is a rather strong 
in-plane uniaxial magnetic anisotropy of epitaxial 
(Ga,Mn)As layers grown on GaAs substrates of (001) ori- 
entation. Owing to the presence of the twofold symmetry 
axes [100] and [010], the in-plane zinc-blende directions 
[110] and [TlO] are expected to be equivalent. Yet, as im- 
plied by the character of magnetic anisotropy, the sym- 
metry is lowered from Did to Ci v , possibly due to the 
growth-induced lack of symmetry between the bottom 
and the top of the layer, 9,10 ' 16 which can be phenomeno- 
logically described by introducing a shear strain. 10,17 ' 18 

Since these symmetry considerations are limited to 



(001) layers, investigation of layers grown on substrates of 
other orientations may not only allow to compare exper- 
imental observations with predictions of the p — d Zener 
model in a more general situation, but also provide in- 
formation from which conclusions on the nature of the 
additional anisotropy can be drawn. 

In this paper we present results of studies on 
Gai^^Mn^As layers grown by low-temperature 
molecular beam epitaxy (MBE) on GaAs substrates 
with the (113)A orientation. Previously, magnetic 
anisotropy in such films was probed at low tem- 
peratures by magnetoresistance 19-21 , scanning Hall 
probe microscopy 22 and ferromagnetic resonance 
measurements. 21 ' 23,24 Experimental techniques em- 
ployed here include superconducting quantum inter- 
ference device (SQUID) magnetometry, ferromagnetic 
resonance (FMR), and polar magnetooptical Kerr effect 
(PMOKE). Our measurements are carried out over a 
wide temperature and magnetic field range. We find 
that magnetic anisotropy can be consistently described 
taking into account three contributions: a uniaxial 
anisotropy with the hard axis tilted from [113] toward 
the [001] direction, an in-plane uniaxial anisotropy with 
the easy axis along the [TlO] direction, and a cubic 
anisotropy with easy (100) directions. The general form 
of anisotropy is, therefore, similar to the case of (001) 
films, but the direction of the hard axis is found to be 
neither along [001] nor perpendicular to the film plane 
in the (113) case. The accumulated experimental results 
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allow us to determine how the three relevant magnetic 
anisotropy constants K as well as the tilt angle depend 
on the temperature. We find that the magnitudes of 
energies corresponding to the competing cubic and 
uniaxial anisotropies in the (001) plane depend, as 
could be expected, as the fourth and second power 
of spontaneous magnetization M(T), respectively. In 
contrast, a complex dependence on M(T) is observed in 
the case of the energy characterizing the out-of-plane 
uniaxial anisotropy. We assign this behavior to the 
spin-splitting-induced and, hence, temperature depen- 
dent redistribution of holes between the valence band 
subbands that are characterized by different directions 
of the angular momentum and, hence, of the easy axes. 

In the theoretical part, we present a theory of magnetic 
anisotropy in epitaxially strained layers of (Ga,Mn)As 
and related systems within the p — d Zener model. Our 
approach generalizes earlier theories developed for (001) 
films 4 ' 5 ' 17,18 ' 25 by allowing for an arbitrary crystallo- 
graphic orientation of the substrate. Similarly to pre- 
vious studies, 10 ' 17 in order to explain the experimental 
findings, we introduce an additional shear strain whose 
three components constitute adjustable parameters. We 
also take into account the Hamiltonian terms linear in 
k and find that they give a minor contribution to the 
magnitude of magnetic anisotropy constants. 

II. SAMPLES AND EXPERIMENT 

We study a 50 nm thick Gai-^Mn^As layer which has 
been grown on a (113) A GaAs substrate (see Fig. 1) by 
low temperature MBE. 26 The total Mn concentration of 
x = 6.4% has been determined by means of secondary ion 
mass spectrometry, however only more than twice lower 
value of an effective Mn concentration x e g can be inferred 
from low temperature experimental saturation magneti- 
zation, Af cxp . This reduction of x is primarily caused 
by a presence of Mn interstitials. These point defects 




FIG. 1: (Color online) Crystallographic directions for a GaAs 
substrate of (113) orientation. 



act as double donors and form strongly coupled spin sin- 
glet pairs with neighbor substitutional Mn cations. 6,27 ' 28 
These pairs neither participate in the ferromagnetic or- 
der, nor they contribute to M. Thus, the effective con- 
centration of Mn ions which generates M oxp gets reduced 
to x e s — x — 2xi, where x\Nq is the concentration of 
the Mn interstitials and Nq is the cation concentration. 
However, the experimentally measured M cxp is further 
reduced by holes magnetization, Mh, which is oppositely 
oriented to magnetization of Mn spins, Mm,,, and so My^ n 
= M exp + \Mh\ should be used to calculate x c g, with 
Mh being computed in the framework of the the mean- 
field p — d Zener model. 4,5 We perform these calculations 
in a self-consistent way taking the hole concentration as 
p = Nq(x — 3x\) = No(3x c s — x)/2, that is neglecting 
other charge compensating defects. 

The open air post growth annealing at temperatures 
below or comparable to the growth temperature 29 ' 30 is a 
frequently used procedure for improving material param- 
eters of (Ga,Mn) As, since the corresponding out-diffusion 
and passivation of Mn interstitials 31 increases x e g, p, and 
eventually the Curie temperature Tq. Therefore in order 
to widen the parameter space employed here to study 
the magnetic anisotropy in this compound we investigate 
both the as-grown material (sample SI) and the samples 
annealed at 200°C for 1.5 hour (sample S2) and 5 hours 
(sample S3). Taking the determined values of M exp data 
we end up with x e s = 2.7, 3.1, and 3.3% and p = 2.0, 
3.3, and 3.8 x 10 20 cm~ 3 for which calculated values of 
Tq = 46, 73, and 85 K compares favorably with the ex- 
perimentally established values of 65, 77, and 79 K, for 
samples SI, S2, and S3, respectively. 

Magnetic properties referred to above and described 
further on have been obtained by employing a Quan- 
tum Design MPMS XL-5 magnetometer. A special de- 
magnetization procedure has been employed to minimize 
the influence of parasite fields on zero-field measure- 
ments. The temperature dependence of remnant magne- 
tization, TRM, serves to obtain an overview of magnetic 
anisotropy as well as to determine Tq (Sec. Ill A). After 
cooling the sample across Tc down to 5 K in the exter- 
nal magnetic field of 0.1 T, the field is removed, allowing 
magnetization to assume the direction along the closest 
easy axis. The magnitude of the magnetization compo- 
nent along the magnet axis, TMRj, is then measured 
while heating, where i indicates one of the three mutu- 
ally orthogonal directions of the magnetizing field, corre- 
sponding to the surface normal ri\ = [113] and the two 
edges n 2 = [332], n.3 — [110], as depicted in Fig. 1. Since, 
except to the immediate vicinity of the spin reorienta- 
tion transition, magnetization of (Ga,Mn)As films tend 
to align in a single domain state, the measurements per- 
formed for the three orthogonal axes provide the temper- 
ature dependence of the magnetization magnitude and 
direction. 

To study magnetic anisotropy in a greater detail, mag- 
netic hysteresis loops M^(if) have been recorded in ex- 
ternal magnetic field in the range of ±0.5 T along the 
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three directions i. The measurements have been car- 
ried out at various temperatures, and the parameters 
of the anisotropy model (Sec. IIIB) have been fitted 
to reproduce the magnetization data. To cross check 
magnetic anisotropy constants obtained from SQUID 
studies, FMR measurements have been performed at 
uj/2tt = 9.3 GHz and T = 10 K. We have performed 
angle dependent measurements of the resonance field in 
the four different crystallographic planes (110), (332), 
(113) and ^^(3 - VTT, 3 + VTT, -2). As discussed in 
Sec. IIIC, the FMR data are in a good agreement with 
the anisotropy model, employing parameters determined 
from the SQUID measurements. 



III. EXPERIMENTAL RESULTS 
A. Overview of magnetic anisotropy 

The TRM studies of all three samples are summarized 
in Fig. 2. We immediately see that the TRMp- 1Q j com- 
ponent of TRM is the strongest for all of the samples 
and that at elevated temperatures its magnitude is nearly 
equal to the saturation magnetization M(T), established 
by the measurement in fioH = 0.1 T. Since the magni- 
tude of the other two magnetization components is van- 
ishingly small, we find that in this temperature range 
the in-plane uniaxial anisotropy with the easy axis along 
[TlO] direction dominates. This perfectly uniaxial behav- 
ior at T — > Tq allows us to use TRMj- 10 j to precisely 
determine Tq in the studied samples (already given in 
the previous section). On the other hand, below a cer- 
tain temperature T* (marked by arrow for every sam- 
ple in Fig. 2) TRM^q, gets visibly smaller than M(T), 
and the other in-plane TRM component, TRMj 33 2j, ac- 
quires sizable values, followed at still lower temperatures 
by the out-of-plane component TRMr 113 i. This clearly 
indicates a departure of the easy direction from the [TlO] 
direction below these characteristic temperatures. Such 
a scheme turns out to be fully equivalent to the general 
pattern of magnetic anisotropy in (001) (Ga,Mn)As un- 
der compressive strain. 7 ' 8 ' 10,11 ' 16 ' 32 ' 33 In such films uni- 
axial anisotropy between [110] and [110] directions, dom- 
inating at elevated temperatures, gives way at low T to 
biaxial anisotropy with in-plane (100) easy axes. This 
spin reorientation transition takes place at a tempera- 
ture, at which uniaxial and biaxial anisotropy constants 
equilibrate, 11 and is corroborated numerically in our sam- 
ples from analysis of the magnetization processes pre- 
sented in Sec. IIIB. 

In an analogy to (001) (Ga,Mn)As, let us assume for 
a moment that M of a (113) sample remains (without 
a magnetic field) in (001) plane. Then, a similar de- 
scription in terms of two in-plane anisotropies (one biax- 
ial and one uniaxial) is possible. Furthermore, assuming 
that the uniaxial anisotropy constant is proportional to 
M(T) 2 , the biaxial anisotropy constant is proportional 
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FIG. 2: (Color online) Temperature dependence of rem- 
nant magnetization components in all three samples (points). 
Black lines: M(H) at 0.1 T along [TlO] - the in-plane uniaxial 
easy axis. Solid lines of lighter shades: the same magneti- 
zation components calculated according to a model of only 
two: uniaxial and biaxial magnetic anisotropies operating at 
(001) plane and undergoing a spin reorientation transition at 
temperatures T* (marked by arrows), as in (001) (Ga,Mn)As. 



to M(T) 4 (Ref. 11), and that both equilibrate at T* we 
are able to model qualitatively temperature induced ro- 
tation of magnetization in the sample and calculate all 
three components of magnetization that would be mea- 
sured by SQUID. The thick solid lines in Fig. 2 show 
the results for sample SI and we find them reproducing 
the experimental findings reasonably well. Therefore we 
identify T* as the temperature at which the spin reorien- 
tation transition from a biaxial anisotropy along (100) to 
uniaxial one along [110] takes place in this system. On 
the other hand, the discrepancies seen in Fig. 2 indicate 
that a more elaborated model is needed. In particular, 
we can infer from the low temperature TRM data that 
the orientation of M at 5 K moves actually away from 
(100) on annealing. The angle between M and (100) 
is increasing from 9, through 19 to 26 deg for samples 
SI, S2 and S3 respectively. At the same time the an- 
gle between M and (113) plane is dropping from 13 to 
7 deg. This indicates that the plane in which both easy 
orientations of M reside at low T is tilting away from 
(001) towards (113) plane. This observation is fully con- 
firmed from the comprehensive analysis of the magnetic 
anisotropy presented in the next section. 

We remark here that the origin of the symmetry break- 
ing between [110] and [lTO] in (001) (Ga,Mn)As is still 
unknown and it is very stimulating to see favoring the 
in-plane [110] also in layers of different surface recon- 
struction than (001) GaAs. 
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B. Experimental determination of anisotropy 
constants 

In order to build up a more complete anisotropy de- 
scription we analyze full magnetization curves M(H). 
It was shown by Limmer et al. (Ref. 20,21) that an 
accurate description the magnetic anisotropy in (113)A 
(Ga,Mn)As requires at least four components: a cubic 
magnetic anisotropy with respect to the (001) axes, uni- 
axial in-plane anisotropy along the [110] direction, and 
two uniaxial out-of-plane anisotropies along the [113] and 
[001] directions. The first two anisotropy components are 
commonly observed in (OOl)-oriented (Ga,Mn)As sam- 
ples, and as shown in the previous section they are suf- 
ficient to provide even semi-quantitative description in 
(113) case. The other two arise from the epitaxial strain 
and demagnetizing effect, both of which depend on the 
orientation of the substrate. In our approach we combine 
the two out-of-plane magnetic anisotropy contributions 
into a single one, with its hard axis oriented between the 
[001] and [113] directions. Accordingly, we write the free 
energy in the form, 



F = -fj, H ■ M + Kc{w 2 x w 2 y +w 2 y w 2 z +w 2 z wl) 
+ K uTw sin2 6 sin2 * 

+ K u \ (cos ©a cos — sin ©a sin© cos $) 2 . (1) 



Here, Kn, K„ 



and K,n are the lowest order cubic, in- 



plane uniaxial and out-of-plane uniaxial anisotropy en- 
ergies, respectively; 0a describes the angle between the 
K u \ hard axis and [113] direction; w x , w y , and w z de- 
note direction cosines of the magnetization vector with 
respect to the main crystallographic directions (100); 
is the angle between M and the [113] direction, and $ 
is the angle between the projection of the M onto the 
sample plane and the [332] direction. 

By numerical minimizing of the free energy with re- 
spect to and <E> we are able to trace the rotation of 
M, starting from the given orientation, while sweeping 
or rotating external magnetic field. Adjusting the ob- 
tained "trace" to the experimental data we get the val- 
ues of the four parameters of the model. We perform 
this procedure numerically for every sample for all three 
orientations and for all temperatures the Mi(H) curves 
have been recorded. Figure 3 shows an example of the 
measured and fitted Mi(H) for the sample SI at two dif- 
ferent temperatures. 

Temperature dependence of the three magnetic 
anisotropy constants and the angle 0a are presented in 
Fig. 4a and b, respectively. All Ki's monotonically de- 
crease with temperature, and, like in (001) (Ga,Mn)As, 
the cubic anisotropy constant Kc (with (100) easy axes, 
see Fig. 4c) and in-plane uniaxial constant K u j 1Q (with 
[110] easy axis, see Fig. 4d) are proportional to M 4 and 
M , respectively, so confirming the validity of the sin- 
gle domain approach used to analyze the observed mag- 
netization rotations. The Kq{T) and K u j 1Q (T) data 
point to the presence of the spin reorientation transi- 





30 




20 - 


E 






10 - 


C 




Q 


- 


CD 




.y 




"(D 


-10 - 


cn 




CD 






-20 - 




-30 - 




30 r 



(a) T= 10K 




□ 


□ 


> 




exp. fitting 
results results 



[110] 
[332] 
[113] 



0.0 0.1 0.2 
M H (T) 



0.4 



FIG. 3: (Color online) Examples of magnetization curves for 
the as-grown sample measured along three, mutually orthogo- 
nal, major sample directions [TlO], [332], and [113]: (a) below 
spin reorientation transition at 10 K and (b) above, at 38 K. 
Symbols indicate measurement points; lines represent the best 
fit of the model described by Eq. (1). 



tion in (001) plane. This already inferred from TRM 
data magnetic easy axis changeover must take place as 
Kc(T) and K u j w (T) swap their intensities in our sam- 
ples. The relevant temperatures are marked in Fig. 4a 
by arrows. Importantly, we find these temperatures to 
agree within 2-3 K with those indicated in Fig. 2, what 
strongly underlines the correctness of the approach we 
employ here to describe the magnetic anisotropy in our 
samples. We note that the SRT shifts to lower tempera- 
tures on going from sample SI to S3, since on annealing 
the in-plane uniaxial anisotropy gets strongly enhanced 
relative to the cubic one (compare Figs. 4c and d). 

In contrast, K u \ shows a more complex dependence 
on M 2 , see Fig. 4e. A proportionality of the out-of-plane 
anisotropy constant to M 2 is seen only at low M, that 
is at high T. On lowering temperature K u \ departures 
from this trend, and the effect is the strongest for the 
sample SI. We connect this behavior with the proximity 
of this system to another spin reorientation transition, 
the transition from the hard to easy out-of-plane axis 
of the K u \ uniaxial magnetic anisotropy. Such a SRT 
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FIG. 4: (Color online) (a) and (b) points: temperature dependence of Kc, K u j 10 and K u \ (a) and angle Qa (b) obtained from 
numerical fitting of Eq. (1) to experimental magnetization curves for all three samples considered in this study. Solid, dashed 
and dotted arrows in (a) indicate the spin reorientation temperature in samples SI, S2 and S3, respectively, (c), (d) and (e): 
Kc, K u j 10 and K u \ dependence on M 4 , M 2 and M 2 , respectively. In all panels, the various lines are guides for eye only. 



takes place in compressively strained (001) (Ga,Mn)As 
on lowering T, and was already observed in samples with 
moderate or high x but rather low hole density 9 . The 
effect depends on the ratio of valence band spin split- 
ting to the Fermi energy. Therefore the SI sample, the 
one with the lowest p is expected to show the strongest 
deviations from the expected functional form. Then on 
annealing, along with the increase of p, we expect the 
so called in-plane magnetic anisotropy (for the compres- 
sively strained layers) to become more robust [less de- 
pendent on the magnitude of the valence band splitting, 
that is on M(T)], as experimentally observed. 

Finally, we comment on Qa, the parameter that can 
serve as a measure of the angle between an 'easy plane' 
with respect to K u \ (perpendicular) hard axis and the 
sample face. As indicated in Fig. 4b that angle remains 
nearly constant at elevated temperatures and shows a 
weak, but noticeable turn towards [001] below tempera- 
tures which can be associated with K u j w <^> Kc SRT. 
This behavior again indicates the departure of easy direc- 
tion of M from [110] direction in the (113) plane. How- 
ever, the maximum determined value of Qa — 15° indi- 



cates, that the rotation of M actually neither takes place 
in the (001) plane, nor is it directed exactly towards (100) 
directions. M rather follows a complex route in between 
(001) and (113) planes, a conclusion that is a numerical 
confirmation of the results of the simple analysis of the 
TRM data presented in the previous section. 



C. Ferromagnetic resonance 

A tool widely used to study magnetic anisotropy is fer- 
romagnetic resonance spectroscopy. Magnetic anisotropy 
in thin (Ga,Mn)As films on (113)A GaAs was recently 
studied by Bihler 23 and Limmer 20,21 . In a ferromagnetic 
resonance experiment the magnetization vector M of 
the sample precesses around its equilibrium direction in 
given external magnetic field H with Larmor frequency 
lol- The resonant condition at a fixed frequency of mi- 
crowaves uj is given by, 
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FIG. 5: (Color online) Angular dependence of the ferromag- 
netic resonance fields for external magnetic field rotating in 
four different crystallographic planes: (113), (110), (332) and 
j-j=(3 — vTT, 3 + \/TT, — 2). Points show measured reso- 
nance field values. Lines show resonant field values calculated 
using anisotropy energy described by Eq. (1) and magnetic 
anisotropy constants obtained from hysteresis loops. 

Here, 7 = g^s^ 1 is the gyromagnetic ratio, g is the g- 
factor, hb the Bohr magneton, and h is the Planck con- 
stant. Resonance field is obtained by evaluating Eq. (2) 
at the equilibrium position of M (dF/d<d = and 
dF/d$ = 0). 

In Fig. 5 the dependence of the measured resonant 
fields on the orientation of the applied magnetic field 
is shown for the sample S2, along with the results of 
a calculation made according to Eq. (2) with the mag- 
netic anisotropy parameters obtained from SQUID mag- 
netization curves. The agreement between the calcula- 
tion and the measured data is very good, indicating that 
Eq. (1) captures main features of magnetic anisotropy 
and that the numerical procedure employed to extract 
the anisotropy constants is correct. 



D. Magnetization reversal 

We end the experimental part evidencing an interesting 
mechanism of the reversal of the out-of-plane magnetiza- 
tion component by an in-plane magnetic field. 

Below the spin reorientation transition (about 20-30 K) 
the magnetization easy axes are moving close to the (100) 
directions, that is they are tilted up from the sample face, 
the (113) plane, and so M is acquiring a sizable non-zero 
component Mr 113 ]. These two axes define a plane lying 
between (113) and (001) planes, which shares the com- 
mon [TlO] direction with those two planes. Therefore, any 
sweep of an external field, except that along the [110] di- 
rection, will result in a magnetization rotation across the 
[TlO] line from one half of that plane (say that one 'above' 
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FIG. 6: (Color online) The out-of-plane magnetization com- 
ponent dependence on the in-plane magnetic field (swept 
along [332]). The measurement is performed on the as-grown 
sample SI at T = 15 K. The cartoons inserted in the figure 
illustrate the process. 



sample face) to the other one (say 'below') resulting in 
Mr 113 i reversal. 

Such a process is illustrated in Fig. 6 in the most inter- 
esting case, when the field is swept in the sample plane, 
along [332]. We record the Mr 113 i magnetization compo- 
nent using the polar MOKE technique and a clear change 
of sign of the signal evidences the reversal of M by the 
application of an in-plane magnetic field. The cartoons 
inserted in this figure visualize the mechanism of this re- 
versal. 



IV. THEORY OF MAGNETIC ANISOTROPY 

A. k ■ p Hamiltonian 

The current theory describing the properties of the 
(Ga,Mn) As ferromagnetic semiconductor is the p-d Zener 
model. 4 In this model, the thermodynamic properties are 
determined by the valence band carriers contribution to 
the free energy of the system, which is calculated taking 
the spin-orbit interaction into account within the k ■ p 
theory 4-6 ' 17 or tight binding model 34 with the p-d ex- 
change interaction between the carriers and the local- 
ized Mn spins considered within the virtual-crystal and 
molecular-field approximations. Within this approach, 
magnetic anisotropy depends on the strain tensor com- 
ponents. 

The 6-band Luttinger k ■ p Hamiltonian of a valence 
band electron in a zinc-blende semiconductor is a block 
matrix (cf. Ref. 35): 

_ fH vv H*>\ 

''•6x6 — I Jjsv JJ SS ] ' v / 
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where; 



m\2 n 



72 



J x — J I k„ 



c.p. 



- 273 [{J x , Jy}{k x , k y } + c.p.}} 

h 2 

-— [-il2{U xx k 2 x + c.p.) 

- 6-f3(U xy {k x ,k y } + c.p.)} 

2m 



(4) 

(5) 
(6) 



(we use the notation of Ref. 35). Our basis is related to 
that of Ref. 5 as follows: u\ — — ||, |), u 2 = — i ||, |), 

U3 = I § , — ^> , Ui = i|f,-f), "5 = M 6 = 

i ||, — |), i.e. we use the standard basis of angular mo- 
mentum eigenvectors (notice the change of sign in ||, 5) 
and || , — |) with respect to Ref. 35 that accounts for the 
difference in the sign of H vs , H sv ). In this basis the p — d 
exchange Hamiltonian is: 



U pd = B G 



2(J-w) -6(U-w) 
-6(T-w) -(a-w) 



(7) 



where w 
Ref. 5, 



M/M and Bq is given by equation 2 of 



B G = A F pM/{6gti B ), 



while the strain Hamiltonian is: 



nr = -b 



(Jx - 7^J 2 )txx + C.p. 



H v e s 



d 

"71 

= -3b(U xx e xx +c.p.) 



V3d(2U xy e xy + c.p.) 



0. 



(8) 

(9) 
(10) 

(11) 
(12) 
(13) 



Since the strain tensor for the (113) substrate orien- 
tation features non-zero non-diagonal components, it is 
necessary to include in the k p Hamiltonian the so called 
fc-linear terms, i.e., terms linear in k and e coming via 
second-order perturbation (Ref. 36, §15) from the terms 
in the 8x8 Kane Hamiltonian 37 that mix the conduction 
and the valence bands. The corresponding Hamiltonian 



is 



Hke = C4 



(1 



J ■ (f 



f(l-i)^ (14 ) 



T ■ (f (1 — ry) er • ip 



where 77 = Ao/(E g + Aq) and the components of the 
vector cp are ip z = e zx k x — e zy k y (c.p.). The numerical 
value given in Ref. 37 is Cs/h = 8 x 10 5 ms _1 , hence for 
C 4 = -63/(277) we obtain C A /h = -2.2 x 10 6 m/s. 



B. Strain tensor 

Determining the components of the strain tensor for an 
unrelaxed epitaxial layer grown on a lattice mismatched 
substrate can be considered a classical topic. The two 
possible approaches to this problem are (i) to solve a 
system of linear equations for the strain and stress com- 
ponents assuming that some components of those ten- 
sors vanish (this is our approach) or (ii) to determine the 
strain of the layer by minimizing the elastic energy (this 
is the approach formulated in Ref. 38). Our approach in- 
volves a transformation of the coordinate system that is 
feasible in general only using a computer algebra system. 
Using one we arrive to the form of the symmetric strain 
tensor that is in a perfect agreement with that of Ref. 38. 

The result for a (lln) oriented substrate is, 



/ 



A + B + C 



'B + C -A -D/2\ 

-A B + C -D/2 , (15) 
-D/2 -D/2 B-2C) 



where 



A = 3 [n 2 (c n - ci 2 ) - (n 2 - 2)c 44 ] (cu + 2c 12 ), (16) 

B = 2 [(n 4 + 4n 2 + l)(c n - c 12 ) + 6c 44 ] c 44 , (17) 

C = (n 2 -l)(n 2 + 2)( Cll + 2c 12 )c 44 , (18) 

D = 3n [(n 2 + l)(c n - c 12 ) + 2c 44 ] (c n + 2c 12 ), (19) 



and 



/ = -Aa/a = (a - a)/c 



(20) 



is the relative lattice constant misfit. Here, the compo- 
nents of the epitaxial strain tensor are given with respect 
to the coordinates (x, y, z) associated with the crystallo- 
graphic axes, x = [100], y — [010], and z = [001] (i.e. the 
quantities given above are the strain components that 
enter the k ■ p Hamiltonian) . 

Using the values Cn = 119. GPa, ci 2 = 53.8 GPa, and 
c 44 = 59.5 GPa (Ref. 39, p. 105) we obtain 



-cpi 



0.9488 -0.0512 -0.3478\ 
-0.0512 0.9488 -0.3478 
k -0.3478 -0.3478 -0.6260/ 



(21) 



For the purpose of determining the strain components 
from X-ray diffraction data, the components of the strain 
tensor in the coordinate system associated with the epi- 
taxial film are needed. We take as the coordinate system: 
x' = [3,3,-2], y' = [-1,1,0], z' = [1,1,3]. The relative 
difference of the lattice constants along the [113] direc- 
tion between that layer and the substrate is 
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Ad/d = [27(cn-c 12 ) + 67c 44 ](c 11 +2c 12 ) = 

7 9(c 11 -c 12 )(c 11 +2c 12 ) + (101cii-34c 12 +4c 4 4)c44 ' ' V ' 

For the sake of completeness we notice that there is also a shear strain component 

12>/2[(cii - ci 2 ) - 2c 44 ](cii + 2ci 2 ) 



9(cu - ci 2 )(cn + 2c i2 ) + (lOlcn - 34c i2 + 4c 44 )c 44 



Ao/a = -0.2746 Ao/o. (23) 



Following Ref. 10, to account for the mechanism which 
generates the in-plane uniaxial anisotropy in (001) sam- 
ples, we incorporate in the p — d Zener model an addi- 
tional Hamiltonian term corresponding to shear strain 



^epi ~t~ ^ • 



(24) 



In case of a (OOl)-oriented substrate the additional strain 
e' has a non-zero xy component, e' xy . The correspond- 
ing anisotropy is of the form K xy w x w y (as in Ref. 40), 
hence it is a difference of uniaxial anisotropies on the 
[110] and [TlO] directions, and the anisotropy field is 
H u = 2K xy /([i M) (this is the field required to align 
magnetization along the hard axis, e.g. [110]; only H u /2 
is required to align magnetization along the z direction). 
In the case of a (113)-oriented substrate the additional 
strain may have more non-zero components. We assume 
that the mirror symmetry with respect to the (110) plane 
is preserved, hence e' xz = e' yz . 



C. Numerical procedure 

The numerical procedure serving to determine the 
magnetic anisotropy from the Hamiltonian matrix is de- 
scribed in Ref. 5. Let us note that including the fc-linear 
terms in the Hamiltonian leads to a tenfold increase of the 
processing time, although in specific cases it is possible 
to generate a symbolic expression for the characteristic 
polynomial of the 6x6 Hamiltonian matrix. Moreover, 
since numerical interpolation of the dependence of the 
hole concentration on the Fermi energy may lead to un- 
controllable inaccuracies, an alternative procedure that 
avoids those inaccuracies is to directly integrate the en- 
ergy of the carriers in the momentum space. However, 
the integration has to be done separately for each hole 
concentration (this is an advantage if a single hole con- 
centration is specified). Moreover, one still needs to solve 
the inverse eigenvalue problem to find the discontinuities 
of the integrand. 

In a numerical calculation, it is possible to determine 
the full magnetic anisotropy by computing the free energy 
of the carriers for a number of directions of magnetiza- 
tion. In our case we choose a grid of directions that is 
rectangular in the spherical coordinates (w x = sin 9 cos 6, 
w v = sin 9 sin 6, w z = cos 9), i.e. 9 = 9; and 6 = 6a, 



where cos^i, i — 1, 2, . . . , Ng are the nodes of a Gaussian 
quadrature and 6j — 2nj/N ( f > , j = 0,1,..., N<p — 1 are 
equally spaced. Then, following the method used in the 
software package SHTOOLS 41 (routine SHExpandGLQ) we 
expand the magnetic anisotropy (free energy) into a sum 
of low-order spherical harmonics. Since the free energy is 
even, choosing even allows to restrict the grid to a half 
of the sphere. We use the standard quantum mechanics 
(orthonormalized) spherical harmonics Yi m (9, 6), and de- 
note the coefficients of this expansion r; m , m = 0, 1, . . . , I, 
and si mi m = 1,2,... I, 



m—1 



(25) 



where the outer sum is over I = 0, 2, . . . , N$ — 1. The 
scalar product in this representation is diagonal, with a 
weight of 1 for r; and 1/2 for r/ m , s; TO , m > 0. 



D. Magnetic anisotropy 

There are a few sources of magnetic anisotropy in epi- 
taxial Gai^Mn^As: the cubic anisotropy of the va- 
lence band, epitaxial strain, the additional off-diagonal 
strain e' , and the shape anisotropy caused by the demag- 
netization effect. Since e' is unknown, it is inevitable 
to parametrize the anisotropy in a manner that sepa- 
rates the components affected by e' from what is pre- 
dictable. If we measure (9, <fi) in the spherical harmonic 
representation of the carriers' free energy with respect 
to the crystallographic axes, the above spherical har- 
monic representation allows one to separate the compo- 
nents brought about by a non-zero value of e' from the 
remaining sources of magnetic anisotropy. Indeed, as far 
as I — 2 is concerned, e' affects primarily only K xy , K yz , 



and K x 



where K xy = 



--S22 



K yz = -, 



:s 2 i, and 



K xz = 
r 22 . As r 22 



|f^2i ■ The remaining components are r 2 o and 
due to the mirror symmetry, they can 
be collected into one term K u00 i(wl — 3), with if„ooi = 

-r 2 o- Thus, we describe the I = 2 anisotropy (with- 
out the demagnetization contribution) by K xy , K xz = 
K yz , and K u qq\. Finally, the cubic anisotropy corre- 



sponds to (r 40 , r 44 ) 



2y^F 
15 



K C {1, y/W?)- 
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We have to relate now the components of the spher- 
ical harmonic representation, r; m and S; m) to the 
experimentally determined magnetic anisotropy con- 
stants Kc, K u x, K u yi , and ®a, as specified in Eq. 1 
and presented in Fig. 4. Since the demagnetization 
effect adds to (7*20, ^21, S21, ^22, S22) the contribution 
^|pX d (4V2, -3V3, -3Vi, 0, v%), with iCj = a*oM 2 /2, 
the constants are related to those of Eq. (1) as follows, 



Kuooi — 



l + 3cos2e^ I 8 

K ul - -A„ Tl0 - — Ad ' i 26 ) 

1- cos 2& A 2 

K ul - A uTl0 - — A d , K 2 ') 



K xz = - 



y2sin2e' 4 



K 



»i 



6 

11" 



-K d 



(28) 



where Q' A — 6^4 — arccos(3/vTT). 

We carry out numerical calculations with band struc- 
ture parameters and deformation potentials specified 
previously. 5 ' 10 We include hole-hole exchange interac- 
tions via the Landau parameter of the susceptibility en- 
hancement, Ap — 1.2 (Ref. 5). This parameter, assumed 
here to be independent of the hole density and strain, 
enters into the relation between M and Be, but also 
divides the anisotropy constants. More specifically, we 
make the calculation with Bq enhanced by the factor Ap, 
and divide the resulting anisotropy constants by Ap -1 , 
where n is the power of magnetization M to which a given 
anisotropy constant is proportional. The result is propor- 
tional to Ap. We have n — 2 for the uniaxial anisotropics 
and n = 4 for the lowest order cubic anisotropy (the pro- 
portionality holds for Bg smaller than a few meV). We 
note that the cubic anisotropy field shown in Fig. 9 of 
Ref. 5 was divided by Ap rather than Ap. 

To evaluate the effect of the fc-linear terms, we use 
a non-zero value of C4 and calculate the difference of 



the resulting anisotropy with respect to the C4 = 
case. This difference has only one noticeable component, 
Ar 2 o = r 2 o(Ci % 0) — r 2 o(C4 = 0), which corresponds 
to a uniaxial anisotropy with a [001] axis (or with [100], 
[010] axes for e yz ^ 0, e xz ^ respectively). A plot 
of AX M ooi oc Ar2o is shown in Fig. 7 for e xy = 0.05% 
(as implied by the symmetry, Aif^ooi is second order 
in e xv ). The values are rather small. In fact, assuming 



Aa/a 



0.5%, we have e. 



0.17%, and the mag- 



nitude of Aif „ioo = Aif u oio is below 10 J /m 3 if we con- 
sider the epitaxial strain only. This estimate appears to 
remain valid in case of a general strain of a similar magni- 
tude, although other anisotropy components are affected 
as well and the dependence on strain components is non- 
linear. However, we stress that this estimate depends on 
the value of the parameter C4, which is somewhat un- 
certain and may be different for the ordinary strain and 
e'. Considered this, it is justified to set C4 = in the 
remaining part of this paper. 

Before we proceed to the calculations specific to the 
particular samples, we make a remark that the data orig- 
inally shown in Fig. 6 of Ref. 10 were not correct due to 
a numerical error in the form of the strain Hamiltonian. 
We show corrected results for H u in Fig. 8. The present 
results are in agreement with Fig. 17 of Ref. 17 (remem- 
ber that our model includes the parameter the Landau 
parameter Ap, neglected in Ref. 17). 

As discussed previously, 9 ' 10 owing to sign oscillations 
of anisotropy constants, the direction of magnetization 
can be changed by temperature {Bq) or the hole con- 
centration, particularly in the vicinity of p — 6 x 10 20 
and 1 x 10 20 cm -3 , according to the results displayed 
in Fig. 8 The corresponding in-plane spin reorientation 
transition has indeed been observed by some of us ei- 
ther as a function of temperature 10 or the gate volt- 
age in metal-insulator semiconductor structures 42,43 in 
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FIG. 7: (Color online) The contribution of the fc-linear terms 
to the [001] uniaxial anisotropy for d/h = —2.18 x 10 6 m/s 
and e xy = 0.05%. 



FIG. 8: (Color online) Hole concentration dependence of the 
in-plane uniaxial anisotropy field due to shear strain e xy = 
0.05% for various values of the valence-band spin-splitting. 
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these two hole concentration regions in (Ga,Mn)As, re- 
spectively. 



V. COMPARISON BETWEEN EXPERIMENT 
AND THEORY 

We detailed above a microscopic model of the magnetic 
anisotropy in a DMS. In order to assess the applicability 
of this model to an arbitrarily oriented DMS we com- 
pare its predictions with experimental findings for (113) 
(Ga,Mn)As. 

First, we specify the magnitude of the lattice mismatch 
to establish the components of the strain tensor. Accord- 
ing to Fig. 2 of Ref. 44, for a (113) sample containing 6.4% 
of Mn we expect Ad/d = 5.6 x 10 -3 which, employing 
Eq. 22, translates into Aa/a = 0.323%. We assume this 
value throughout this section. 

Then, using the values of x c g already established in 
Sec. II, we calculate 

M Mn = x cS N SgnB (29) 

and obtain Bq for each of our samples from Eq. 8. It is 
worth repeating here, that the established upon total x 
and M exp values of x e s and p reproduce, within the same 
model, experimental values of Tq remarkably well. This 
boosts our confidence in the accuracy of the material pa- 
rameters used here for the computations of the magnetic 
anisotropy and gives a solid ground to the presented con- 
clusions. 

The calculations are performed as a function of p, em- 
ploying for each sample the corresponding value of Bq: 
-16.4, -18.7, and -19.7 mcV for samples SI, S2, and 
S3, respectively. The results are presented in Fig. 9 as 
curves, while full symbols represent experimentally es- 
tablished values of the anisotropy constants (a square 
for the sample SI, a triangle for S2, and a circle for S3). 
The experimental if„ooi, K xy , K xz , and Kc are obtained 
from K u j 10 , K u i, Kc and Qa using Eqs. 26, 27 and 28. 
We take the T = 5K experimental anisotropy data, as 
this is what is consistent with the T = limit, implicitly 
assumed in Eq. 29. 

We start by discussing the strongest component of the 
magnetic anisotropy, the K u00 i term. The calculated 
curves are presented in Fig. 9a. The calculations have 
been performed without introducing the fictitious shear 
strain e' (i.e. e' = 0). However, since we know that the 
magnitude of K u ooi is negligibly affected by e', the re- 
sults should already match the experimental data, and 
indeed they do. Although the spread of the experimen- 
tal points in Fig. 9a is significantly larger than that of 
the theoretical curves, the correspondence between the 
computed and experimental values is good, and it has 
been achieved without introducing into the model any 
adjustable parameters. This has been only possible by 
including the hole liquid magnetization in the calcula- 
tion of x e s- When Mh is disregarded, the experimental 



values of K u00 i are systematically above the maxima of 
the theoretical curves. 

As already mentioned, in the case of the (001) 
(Ga,Mn)As films an additional low-symmetry term has 
to be introduced into the Hamiltonian in order to re- 
produce the experimentally observed uniaxial in-plane 
magnetic anisotropy. In a general case of an arbitrarily 
oriented substrate, there are three anisotropy constants 
of this kind, K xy , K yz , and K xz . In the case of (113) 
(Ga,Mn)As, the symmetry requires that K xz = K yz , the 
assumption confirmed by the description of the experi- 
mental results. The two relevant anisotropy constants 
K xy and K xz are presented in Figs. 9(b) and (c), and are 
seen to be non-zero even in the absence of a symmetry 
lowering perturbation, e xy — 0. The computed magni- 
tude of K xy for e' xy = yields an acceptable agreement 
with the experimental data. We note, moreover, that an 
exact match is possible when allowing for non-zero values 
of e' xy ^ 0.001% for sample SI and 0.01% for samples S2 
and S3. 

In contrast, according to Fig. 9(c), the theoretical de- 
scription of the experimental values of the K xz anisotropy 
constant requires a quite sizable value of the correspond- 
ing strain component e' xz = —0.1%. Thus, remarkably 
and contrary to the case of (001) (Ga,Mn)As, one barely 
needs any in-plane shear strain to reproduce in-plane uni- 
axial anisotropy, whereas for the out-of-plane component 
two times stronger shear strain is needed comparing to 
the (001) case. This finding should be taken as a strong 
evidence that the in-growth surface reconstruction, and 
a related with it orientational preferences of Mn incorpo- 
ration must play a decisive role in the mechanism leading 
to the lowering of magnetic symmetry. 

Finally, we turn to the case of the cubic anisotropy 
constant Kc, shown in Fig. 9d. Since Kc shows only a 
small sensitivity to e', we present the results of computa- 
tions only for e' = 0. We find that similarly to the (001) 
case, 9 the present theory underestimates the magnitude 
of Kc , particularly in the low hole concentration region, 
where the theoretically expected change of sign of Kc 
is not observed experimentally. The origin of this dis- 
crepancy, and in particular its relation to the symmetry 
lowering perturbation is presently unknown. 

We have also examined theoretically how the particu- 
lar anisotropy constants depend on magnetization M. As 
could be expected, and in a qualitative agreement with 
the experimental finding shown in Fig. 4, the uniaxial 
anisotropy constants -Kuooi, K xy , and K xz (or equiva- 
lently K u \ and K u j w ) are proportional to M 2 , whereas 
Kc to M 4 , except to the hole concentration region in the 
immediate vicinity of the sign change. 



VI. CONCLUSIONS 

We have investigated the magnetic properties of as- 
grown and annealed (Ga,Mn)As layers grown by MBE 
on GaAs substrates of the (113) orientation and provided 



11 



(a) ' fflS3 

x\ •-■•"■•.^ 

/•- S1 

/ (001) easy plane 


e' = 


0.5 

E 

3 


— n 




a 

. ii [001] easy axis 




5 


J 




-0.5 



(b) 



9 ©S3 



e' = 



S1 

[110] easy axis 



J 



[110] easy axis 



E 
—3 

I! 



(c) 

e' =0 



e 1 =-0.1% 

if 

[1 1 V2] easy axis 



[1 1V2] easy axis 



8=0 £ * 

S3 

S2 

S1 







1 



1.5 
1.0 
0.5 
0.0 

y 

-0.5 
-1.0 



~5 



S1 



ft M 

ill <100> easy axis 
1) 

> 



S2 

* i 3 



.. ■'// <1 1 1> easy axis 
V ■■-'/' 
I St 







1 



p(10 cm ) 



p(10 cm ) 



FIG. 9: (Color online) Lines: theoretical dependence of the anisotropy constants (a) -ftTuooi, (b) K xy , (c) K xz , and (d) Kc 
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the most complete to date description of the magnetic 
anisotropy in the whole temperature range up to To 
At higher temperatures the [110] direction is the easy 
magnetization axis before and after annealing. At low 
temperature the spin reorientation transition to a pair of 
easy axes near the [100] and [010] directions takes place 
and to a first approximation, magnetization behavior as 
a function of temperature is similar to that observed in 
(001) (Ga,Mn)As in the absence of an external magnetic 
field. 11 However, the magnetization vector resides in a 
plane close to (001) plane only for a low value of the hole 
concentration. When it increases, the plane rotates along 
[TlO] towards the sample face (113), and the two cubic 
easy directions move towards the [332] direction. 

We have estimated the values of magnetic anisotropy 
constants by fitting our phenomenological model to the 
hysteresis loops measured by SQUID. The comparison to 
results of FMR measurements confirms the correctness of 
this approach. The obtained values of the cubic and uni- 
axial in-plane magnetic anisotropy constants are propor- 
tional to M 4 and M 2 , respectively. Inflections from the 
M 2 dependence of the out-of-plane uniaxial anisotropy 
constant indicate a proximity to another spin reorienta- 



tion transition at which the out-of-plane axis becomes 
easy on lowering temperature. It has been evidenced by 
MOKE that it is possible to reverse the out-of-plane mag- 
netization component by applying an in-plane magnetic 
field. 

For the hole and effective Mn concentrations deter- 
mined from the values of saturation magnetization and 
the total Mn concentration, the p — d Zener model ex- 
plains, with no adjustable parameters, the magnitude 
of the Curie temperature as well as the sign and mag- 
nitude of the uniaxial [001] anisotropy constant K u qqi 
caused by biaxial strain. At the same time, however, 
the predicted values of the cubic anisotropy constant are 
smaller than those found experimentally in the hole con- 
centration range studied here. For the substrate ori- 
entation in question there are two additional non-zero 
second order (I = 2) components K xy and K xz . The 
comparison of their experimental and theoretical values 
points to the presence of an additional shear strain. The 
non-vanishing components of this additional strain are 
e' xz — e' yz w —0.1%, in contrast with (001) samples, for 
which a non-zero value of e' , 0.05% has to be as- 
sumed in order to explain the experimental data. This 
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finding provides a hint that a preferential Mn incorpora- 
tion during the growth process accounts for the mysteri- 
ous lowering of the (Ga,Mn)As symmetry. 
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